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(Dated: January 19, 2013) 

Consider a discrete-time quantum walk on the JV-cycle subject to decoherence both on the coin 
and the position degrees of freedom. By examining the evolution of the density matrix of the system, 
we derive some new conclusions about the asymptotic behavior of the system. When N is odd, the 
density matrix of the system tends, in the long run, to the maximally mixed state, independent of 
the initial state. When N is even, although the behavior of the system is not necessarily asymptot- 
ically stationary, in this case too an explicit formulation is obtained of the asymptotic dynamics of 
the system. Moreover, this approach enables us to specify the limiting behavior of the mutual infor- 
mation, viewed as a measure of quantum entanglement between subsystems (coin and walker). In 
particular, our results provide efficient theoretical confirmation of the findings of previous authors, 
who arrived at their results through extensive numerical simulations. Our results can be attributed 
to an important theorem which, for a generalized random unitary operation, explicitly identifies the 
structure of all of its eigenspaces corresponding to eigenvalues of unit modulus. 

PACS numbers: 03.67.Lx, 05.30.-d, 05.40.-a 



I. INTRODUCTION 

In recent times, quantum walks (QW) have attracted 
extensive attention, mainly for their value as potential 
sources of new algorithms [l]-SJ. However, quantum 
walks are physically and mathematically interesting in 
their own right due to the richness of dynamical and sta- 
tistical properties inherent in the systems 0-0] ■ 

In stride with the theoretical advances, several 
schemes, such as @-[l3|], have been proposed to imple- 
ment quantum walks in realistic media. However, any 
attempt to implement a quantum system in a physical 
channel must take into consideration the critical issue of 
"decoherence". As the aptly contrived neologism sug- 
gests, decoherence connotes disruption of the character- 
istic coherent features of a quantum system, including 
entanglement of subsystems. Subject to decoherence, a 
quantum system tends to behave, in the long-term limit, 
like a classical system. 

Various mathematical models of decoherence in 
discrete-time QWs have been proposed and investigated 
both numerically and analytically [14l - |21| . In general, re- 
gardless of the model adopted, a discrete-time QW, when 
subject to decoherence, tends to revert, in the long-term 
limit, to its classical analogue. However, when prop- 
erly understood and appropriately controlled, it turns 
out that decoherence actually can serve to improve the 
algorithmic properites of a quantum walk [22[ . 

In this paper, we investigate the evolution of a quan- 
tum walk on the A^-cycle under the assumption of 
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decoherence-inducing disturbances on both the coin and 
the position degrees of freedom. As evidenced by the lit- 
erature, much has been written about the case of a purely 
coherent QW on the Af-cycle. For instance, in [23| it is 
shown that if the parity of the cycle- length N is odd, then 
the time-averaged distribution of a coin-governed quan- 
tum walk on the A^-cycle mixes to a uniform distribution. 
On the other hand, according to 0, HI[ , if the parity of 
N is even, then the limiting distribution still exists, but 
might be non-uniform. For an interesting investigation 
of the fluctuation of quantum walks on the A^-cycle, as 
characterized by their temporal standard deviation, the 
reader is referred to [26j | . 

For the QW on the N-cyc\e, the mathematical model 
of decoherence adopted in this paper can be described 
as follows. At every time step of the walk, the option 
persists, with constant probability q, of exercising a pro- 
jective measurement on the basis states of the quantum 
system. As such, three distinct cases emerge: decoher- 
ence might be assumed to apply to 1) the position only, 
2) to the coin only, or 3) jointly to the position and coin. 
Analytic treatments of case 1) can be found in (Tol. [l7j. 
For an analytic treatment of case 2), see [IjIHt]]- Numer- 
ical treatments of all three cases can be found in [HJ [H[ . 
Judging by the literature, analytic treatments of case 3) 
are rare at best. In this paper, our main objective is to 
provide a thorough analytic treatment of case 3). 

By examining the linear structure of the quantum op- 
eration underlying the QW and analyzing the evolution 
of the associated quantum Markov chain, our approach 
yields an explicit formulation of limiting density matrix 
of the quantum system, both for odd and even values of 
the cycle length N. This approach not only allows us 
to derive the limiting probability distribution, but also 
enables us to formulate a clear account of the limiting 
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behavior of the system's quantum entanglement features. 

In what follows, Section 2 is devoted to introducing 
and developing the basic ingredients essential to our in- 
vestigation. In particular, for a generalized random uni- 
tary operation, we specify the linear structure of its key 
eigenspaces. We proceed, in section 3, to examine the 
Markov chain evolution of the density operator of a quan- 
tum walk on the A-cycle subjet to decoherence on both 
the coin and the position degrees of freedom. Our anal- 
ysis includes a consideration of the effect of decoherence 
upon entanglement. In Section 4, we offer some con- 
cluding remarks. Finally, the proofs of the theorems and 
theorem-like assertions are presented in Appendices A 
through E. 

II. GENERALIZED RANDOM UNITARY 
OPERATION 

A well-known model of state transition maps for quan- 
tum systems is provided by random unitary operations 
[29l |3G |. In this section, we introduce a generalization 
of the random unitary operation model. By appending 
to the usual definition of a random unitary operation a 
non-unitary bistochastic component, the resulting model 
is inclusive enough to lend itself to an analysis of the 
long-term tendencies both of decoherence and entangle- 
ment for a quantum walk on the iV-cycle. 

We begin with some preliminaries. Given a Hilbert 
space H of finite dimension m, let 05(H) denote the set 
of all linear operators on H with inner product defined 
by 

(X,Y)=tr(X^Y). (1) 

The corresponding norm, called Frobenius norm or 
Schatten 2- norm, is defined by 

||A|| ee [tr^A)] 1 ^ = yftx^cj. (2) 

This choice of norm on 05(H) will remain in effect 
throughout this paper. 

Let S(H) C 05(H) denote the set of positive operators 
p-.n^U with Tr(p) = 1. The operators p <E D(H) are 
the so-called "density operators". They serve to model, 
as faithfully as do the "state vectors" themselves, the 
possible states of a quantum system whose state vectors 
reside in H. 

By a super- operator <fr on 05(H), we mean a linear 
mapping 3? : 55(H) — > 05(H), with norm defined by 

\\$(X)\\ 

W = SupjfgjB^) . (3) 

Note that dim05(H) = m 2 , where m = dim(H). Thus, 
any super-operator $ on 05(H) can be represented, rel- 
ative to a given basis for 05(H), by an m 2 x m 2 matrix. 
In the sequel, this matrix will be denoted by the sym- 
bol [<&]. In particular, relative to a special basis con- 
sisting of eigenvectors and generalized eigenvectors of 



the shape of the matrix [<&] conforms to a special quasi- 
diagonal lay-out called the Jordan canonical form. The 
details can be found in any one of a number of sources, 
including (3l| . 

Among the set of super-operators, we distinguish a spe- 
cial subset called "quantum operations". By definition, 
to qualify as a quantum operation, the super-operator $ 
must be completely positive, meaning that the extended 
map <fr ® I„ is positive for all n > 1. 

The formalism of quantum operations is versatile 
enough to handle both unitary (closed) and non-unitary 
(open), or a mixture thereof, of discrete transitions of 
state of a quantum system. For a g ood introductory ex- 
position of this subject, see [H, 1 3 3ll 

By Choi's Theorem [34| and [H, [H, HI], any com- 
pletely positive linear operator, including any quantum 
operation $ : 05(H) — > 25(H), can be represented in 
terms of a set A = {Ai \ i = 1, 2, m 2 } of "Kraus oper- 
ators" , as follows: 

*a(X) = Y,A 1 XA\. (4) 

i 

In this expression, which we call the "Choi expansion" of 
3>, the symbol A] denotes Ai (transpose of the complex 
conjugate of Ai). 

In terms of the Choi expansion, the condition of being 
trace-preserving, meaning that Tr($^(A)) = Tr(A) for 
all X G 05(H), is equivalent to the condition: 

E4^=v (5) 

i 

On the other hand, if the Kraus operators of <&yi satisfy 
the dual condition: 

^^4=I ro! (6) 

i 

then is said to be unital. Note that Eq. ([6]) is equiv- 
alent to the simple statement that $^(I m ) = I m . 

A quantum operation which is both unital and trace- 
preserving is called bistochastic. It is a routine matter to 
verify that any convex linear combination of bistochas- 
tic quantum operations is itself a biostochastic quantum 
operation. 

Now suppose the bistochastic quantum operation 3? 
admits a convex decomposition of the form 

k I 

*(p) = E PMP U i + V E A iP A \ > ( ? ) 

i=l 3 = 1 

where the operators Ui are unitary, the operators Aj sat- 
isfy the conditions Eqs. ([5]) and ©, and where, of course, 

q + Yli=i Pi = 1- For lack of a better term, a bistochastic 
quantum operation admitting a convex decomposition 
as in Eq. ([7]) is called a generalized random unitary op- 
eration (henceforth abbreviated GRO). 
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For the purposes of this study, the action of a GRO is 
interpreted as follows. At each time step of the quantum 
walk, the decoherence effect is delivered by the summa- 
tion term in Eq.© involving {Aj} l j =1 , whose coefficient 
q denotes the decoherence rate. Meanwhile, at each time 
step of the walk, the complimentary summation term of 
Eq.© involving Ui imparts a purely unitary transition 
of quantum state whereby each of the Ui is applied to 
the system with corresponding probability pi. 

The GRO model embraces, as special cases, several 
types of quantum operations prevalent in the literature. 
For instance, when q = 0, the GRO model equates to 
the usual model of random unitary operations as in 
[36| . When Ui and Aj are appropriately customized, the 
resulting GRO equates to the model of decoherence, as 
in (l5j . for coin-governed quantum walks. 

Given a bistochastic quantum operation 4>, at least 
one of its eigenvalues A must belong to the unit circle 
(i.e. |A| = 1). In particular, this is so if $ is a GRO 
as in Eq. Q. For an eigenvalue A of <&, let Ker(4> — 
AI) denote the eigenspace of A. Essential to the aims of 
this investigation is a determination of the structure of 
eigenspace Ker($— AI). The following theorem represents 
a significant step in this direction. 

Theorem 1. Let * : S('H) -> <8(H) be the general- 
ized random unitary operation given by 

k I 
1=1 j=l 

and let A, with |A| = 1, be an eigenvalue of 3>. Then 
X 6 Ker(3> — AI) if and only if for each index 1 < i < k 
we have U t X = XXU t and U t XU\ = Y! J=1 AjXA\. 
Proof. See Appendix A. 

Theorem 1 generalizes a result in [3l|. Its proof relies 
on a pattern of reasoning similar to that employed in 
j3lj and the property of the contractiyity of a completely 
positive trace-preserving linear map [37| . 

As shown in literature (e.g. [H, H^]), for a bistochas- 
tic quantum operation, the eigenvalues lying on the unit 
circle determine the evolution of the associated quantum 
Markov process, including the existence or non-existence 
of a long-term stationary state. More precisely, the long- 
term behavior of the quantum Markov process is linked 
intimately to the structure of the eigenspaces of eigenval- 
ues on the unit circle. For an arbitrary quantum opera- 
tion, the problem of determining explicitly all its eigen- 
values on the unit circle remains intractable. However, at 
least for the specific types of quantum operations, the- 
orem 1 provides an efficient means for identifying the 
eigenvalues and the eigenspaces of eigenvalues of absolute 
value 1. The following corollary provides an example of 
this. Nontrivial applications of theorem 1 are relegated 
to the next section. 

Corollary 2. Suppose the bistochastic quantum op- 
eration $ is given by 



^(X) = (l-q)X + q J2^iXAl (8) 

i 

where £\ A t A\ = ql and < q < 1. Then A = I is the 
only eigenvalue on the unit circle. 



III. QUANTUM WALKS ON THE iV-CYCLE 
EXPOSED TO DECOHERENCE BOTH ON COIN 
AND POSITION DEGREES OF FREEDOM 

For a quantum walk on the A-cycle, the position space 
of the walker is the Hilbert space "Hat spanned by an 
orthonormal basis {|x),x € Iin = Z (modiV)}. The coin 
space is the Hilbert space %2 spanned by an orthonormal 
basis {\r}, | I) }. The "state vector space" is % = % at £g) %2 • 
Thus, a general density operator pinH may be expressed 
as 

P= a x ^ yj \xi)(yj\. 

x,y& N i,je{r,l} 

In the sequel, we may use the 2N by 2N matrix [a^yj] 
to represent the density operator p. 

As in [IH , the temporal progression of states of a quan- 
tum walk on the A^-cycle is modeled by repeated itera- 
tions of a quantum operation of the form 

= (1 - q)Uprf +qJ2 E VwiUprfFU. (9) 

xSZjv iG{r,i} 

In Eq. ©, the unitary operator U = S(I<giC) is defined 
in terms of the shift operator S : W — > ~H, which acts on 
basis states \xi) by the formulas S(\xr)) = \{x + l)r) and 
S{\xl)) = \{x — 1)1). As always, I denotes the identity 
operator, acting, in this case, on Hn- Meanwhile, any 
unitary operator C : TL2 —> H-i is eligible to serve as the 
coin operator. In general, C may be expressed by a for- 
mula such as C = un\r) (r\+U2i\l) (r\+ui2\r) (l\+U22\l) 
However, for technical reasons, it will be convenient, in 
the sequel, to assume that the complex coefficients un, 
M12, U21 and U22 ah are non-zero. 

In Eq. (9), the parameter q (0 < q < 1) is called the 
decoherence rate. When q = 0, the QW evolves as a 
purely coherent quantum process. At the other extreme, 
when q = 1, the QW behaves exactly like a classical 
random walk. Since, in this paper, we are interested in 
the case of non-classical quantum walks subject to a non- 
zero level of decoherence, it shall be assumed henceforth 
that < q < 1. Finally, in Eq. the projection 

operator V x i — \xi)(xi\ acts with probability q on the 
Hilbert space spanned by the eigenstate \xi). Thus, at 
each time step of the quantum walk generated by the 
quantum operation exposure to decoherence prevails 
with probability q. 

By comparison with Eq. ((7|) , we see that the definition 
of $ in Eq. §§§ conforms to the definition of a generalized 



4 



random unitary operation (GRO). Thus, <fr falls within 
the purview of Theorem 1 . 

For an eigenvalue A of let E&(\) = Ker(3> — AI) 
denote the eigenspace of A. For 3? as in Eq. ©, the 
following lemma elucidates the relationship between the 
various eigenspaces of 

We digress momentarily to recall what is meant by 
the term "generalized eigenvector" . Suppose J is a Jor- 
dan block in the Jordan canonical form of 4> correspond- 
ing to an eigenvalue A and suppose its dimension is 
dim( J) = m > 1. Let Y\ denote the corresponding eigen- 
vector. Then, by standard linear algebra, there exists a 
sequence Yx, Y 2 , Y m of what are called "generalized 
eigenvectors" characterized by the conditions: 



(&-M)Y k =Y k _x, 



(10) 



where k = l,...,m and where, by definition, Y = 0. 
The generalized eigenvectors belonging to A are linearly 
independent and satisfy the condition (<J? — AI) fe Yfe = 0. 

Lemma 3. Let Ai, A2 be eigenvalues of 3? with |Ai | = 
A2I = 1. Let a be an eigenvalue of 3? with \a\ < 1 and let 
Yx, Yj a denote the generalized eigenvectors belonging 
to a. Then 

1. If Ai ^ A 2 , then E 9 (Xx) 1 £*(A 2 ). 

2. E*(\x) ±Span{Y 1 ,...,Y j J 

Proof. See Appendix B. 

For a QW on the TV-cycle generated by <& as in Eq. 
(O , it turns out that the eigenvalues of <fr on the unit cir- 
cle and the structure of their eigenspaces are completely 
determined by the parity of the cycle length TV. The 
following lemma articulates the details. 

Lemma 4. For a quantum walk on the -/V-cycle gov- 
erned by $ as in Eq. (|9|): 

1. If TV is odd, then 1 is the only eigenvalue on the 
unit circle, and its eigenspace ^(l) = span{I 2 Ar}. 

2. If TV is even, then 1 and -1 are the only eigenval- 
ues of <!> on the unit circle, in which case E&(1) = 
span{l2jv}, while E&( — 1) = span{I±i}, where 
I±i =diag(l, 1,-1,-1..., 1,-1,-1). 

Proof. See Appendix C. 

At this point we have gathered all of the necessary 
ingredients to formulate a description of the long-term 
behavior of a quantum walk on the TV-cycle generated by 
the quantum operation <& as given by Eq.(|9]). 

Theorem 5. Suppose a quantum walk, generated 
by the quantum operation 3? as defined by Eq.([9]), is 
launched on the iV-cycle with initial state p(0) and with 
decoherence rate < q < 1. If TV is odd, then the it- 
erated succession of quantum states p(t) = $ f p(0) con- 
verges to 2]yl2Ar. If N is even, then ||<&'p(0) — jn^2N — 
(— l) t 27v(p(0),I±i)I±i|| converges to zero. In particu- 
lar, if (p(0),I±i) = 0, then p(t) = 4>V(0) converges to 
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Proof. See Appendix D. 

This theorem represents a significant advance com- 
pared to our previous work in [19j. In that paper, which 
dealt with a specific model of quantum walks on the 
TV-cycle, subject to decoherence on the coin degree of 
freedom, an explicit formula was derived for the limit- 
ing probability distribution, not just in the weak time- 
averaged sense, but in the strong point-wise sense (see 
Theorem 3 in [19]). As seen below, this result follows 
as an immediate corollary of Theorem 5 (see Corollary 6 
below). However, unlike Theorem 5 of the present paper, 
Theorem 3 in [l!| fails to specify the limiting structure of 
the density matrices themselves. We speculate that the 
analytical approach employed in this paper might prove 
fruitful also in the context of to specify the limiting 
structure of the corresponding density matrices. 

From this theorem we can derive an immediate corol- 
lary concerning the position probability distribution 
for quantum walks on the TV-cycle. Let P(x, t) = 
Tr (\x)(x\p(t)) denote the probability of finding the 
walker at the position x at time t. 

Corollary 6. For a quantum walk launched on the 
TV-cycle with decoherence rate < q < 1 and driven by 
the quantum operation 3? as in Eq. ©, the following 
assertions hold: 

(i) If TV is odd, then P(x,t) converges to on all 
nodes of the cycle regardless of the initial state p(0). 

(ii) If TV is even and if the quantum walk is launched 
from a definite initial node (i.e., |(p(0),I±i)| = 1) , then 
P(x,t) converges to j| on the supporting nodes of the 
cycle and to on the non-supporting nodes of the cycle. 

(iii) If TV is even and if the parity of the initial node 
of the QW has an equal probability of being odd or even 
(i.e., (p(0),I±i) = 0), then P(x,t) converges to on all 
nodes of the cycle. 

Item s (i) and (ii) of Corollary 6 confirm the predictions 
of [l5l |28| . which are based on numerical simulations. 
However, since the class of quantum walks covered by 
Corollary 6 includes not only those launched in a pure 
state, as considered by [lEIS^, but also those launched 
in a mixed state, Corollary 6 actually is a stronger version 



of the predictions in [15l. |28| 



Our analysis of the long-term evolution of quantum 
walks on the TV-cycle generated by a generalized ran- 
dom unitary operation would not be complete without a 
consideration of the effects of decoherence upon entan- 
glement. We begin by reviewing some preliminaries. 

In what follows, we utilize the concept of von Neumann 
entropy to quantify the mutual information between sub- 
systems (coin and position). Intuitively, the von Neu- 
mann entropy of a quantum system A, denoted S(A), is 
a measure of the uncertainty implied by the multitude 
of potential outcomes as reflected by its density matrix 
p(A). More precisely, S(A) = S(p(A)) = -Tr(plnp). 

For a composite system with two components A 
and B, the joint entropy of their conjunction, de- 
noted by S(A,B), is defined by the formula S(A,B) — 
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— Tr(p AB \np AB ), where p AB is the density matrix of the 
composite quantum system AB. 

A good measure of the level of quantum entanglement 
between the two components A and B is the so-called 
mutual information S(A : B), defined by the formula 
S(A : B) — S{A) + S(B) - S(A, B). 

The following Lemma, due to Watrous [4(| , is essential 
to our reasoning. 

Lemma 7 Let X denote a complex Euclidean space 
and let Pos(A') denote the set of positive semidefinite 
operators defined on X with norm defined by Eq. ([2]l, 
Then, with respect to this norm, the von Neumann en- 
tropy S(p) is continuous at every point p £ Pos(A"). 

For quantum walks on the TV-cycle, the following def- 
initions apply. For the subsystem associated with the 
coin, the time-dependent reduced density operator pc (t) 
is given by pc{t) = tracew(p(i)), where the subscript w 
signifies exclusion or "tracing out" , relative to the over- 
all system density operator p(t), of the walker's degrees 
of freedom. Similarly, for the subsystem associated with 
the walker, the time-dependent reduced density operator 
pw(t) is given by pw(t) — tracec(/o(i)), where the sub- 
script c signifies exclusion or "tracing out" , relative to 
the overall system density operator p(t), of the degrees 
of freedom of the coin. 

The following theorem summarizes our main findings 
in connection with the behavior of quantum entangle- 
ment for quantum walks on the TV-cycle. 

Theorem 8. Suppose a quantum walk is launched on 
the A-cycle with initial state p(0) and with decoherence 
rate < q < 1, and driven by the quantum operation $ 
as in Eq. (|9|). Let pc(t) and pw(t) denote, respectively, 
the time-dependent reduced density operators associated 
with the subsystems of the coin and the walker. Then 
the long-term trend of the mutual information between 
the coin subsystem and the walker subsystem is given by 
lim^oo S (p c (t) : Pw(t)) = 0. 

Proof. See appendix E. 

To put it briefly, according to Theorem 8, exposed to 
any non-zero level of decoherence, the mutual quantum 
information between subsystems tends to zero, signifying 
the total collapse of entanglement between subsystems. 



IV. CONCLUSIONS 

For the bipartite quantum system considered in this 
paper, our analysis has shown that exposure to any 
nonzero level of persistent decoherence causes the sys- 
tem to behave asymptotically like a purely classical sys- 
tem. As noted in [4 11 ] , decoherence on a quantum system 
is manifested through its density matrix by the vanish- 
ing of the off-diagonal elements. Ultimately, the density 
operators should become indistinguishable from diago- 
nal matrices. In the context of a coin-driven quantum 
walk on the iV-cycle, this is precisely what we assert in 
Theorem 5. Note that the off-diagonal elements of the 
density matrix are precisely the elements that represent 



the quantum correlations (also known as entanglement) 
between the coin subsystem and the position subsystem. 
Unsurprisingly, at least for quantum walks, as shown in 
Theorem 8, decoherence turns out to be practically syn- 
onymous with "disentanglement" , so to speak. Indeed, 
as a corollary of our results, we show, when influenced by 
decoherence, the resulting long-term distributions are in- 
distinguishable from those exhibited by classical random 
walks. 

The model of decoherence used in this article is only 
one of several prevalent in the current literature. It would 
be interesting to investigate how quantum entanglement 
and probability distribution respond to other models of 
decoherence (such as amplitude damping, or phase damp- 
ing acting on the coin degree of freedom ), and not just for 
quantum walks on the TV-cycle, but for quantum walks 
over other kinds of topological networks as well. For a 
model of a quantum state transition, the long-term be- 
havior of such associated quantum Markov chains is inti- 
mately linked to the structure of the eigenspaces of eigen- 
values on the unit circle. We speculate that the charac- 
terization of the linear structure of its key eigenspaces in 
Theorem 1, might provide an efficient means for identi- 
fying the eigenspaces of all eigenvalues of absolute value 
1, thereby enabling one to explore the asymptotic evo- 
lution of the associated Markov chains determined by a 
generalized random unitary operation in Eq. ([7]). 
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Appendix A: PROOF OF THEOREM 1 

Proof. If for each index 1 < i < k we have 
UiX = XXUi and U t XU} = Y! j=1 AjXA], then it is easy 
to see that X 6 Ker^— AI). We proceed to justify the as- 
sertion in the opposite direction. Firstly, by [13] , the fol- 
lowing inequality must be true: || Ylj=i AjXAj\\ < \\X\\. 
If X € Ker(5' — AI), then using the unitary invariance 
of the Hilbert-Schmidt norm and the contractivity of a 
completely positive trace-preserving linear map, we get 

k I 

|| A|| = ||AA|| = WY.PiUiXUj + q^jXAlW 

k I 

^pmxujw + qWY^A.xAlw 

1=1 3 =1 
k 

<£ ft ||*|| + g||X|| = |[*ll- (AI) 

i=l 

The inequality (IA1[) actually is an equality. In partic- 
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ular, we have 

l 

11^X^11 = 11X11 (A2) 

For convenience, we set Pk+i = q, = UiXllJ, i = 
l,2,...,k; and v k+1 = Y! 1=1 A 3 XA]. Then Eq. (TATl) 
implies the following equality. 

fe+i fc+i fc+i 

(%2pm,^2piVi) = C^2p i (v l ,v l )^) 2 . (A3) 

z— 1 i—l i—1 

Eq.dAU) leads to 

i<j 

= ^2PiPj[( v hVj) + (Vj,Vi)] 
i<j 

= ^22p iPj Re((v i ,v j )) 

i<j 

< ^22p t p J \(v l ,v :i )\ 

i<j 

<^22p i p j (v h Vi)5(v j ,v j }i. (A4) 

Because of the identity of both right and left hand sides 
of the algebraic expression (jA4|, all inequalities involved 
are in fact are equalities. It can be readily seen that 

Re((«j,Uj)) = \(vi,Vj}\ = (vi,Vi)?{vj,Vj)? ^0 

foxi,j G {l,2,...,fc+l}. (A5) 

According to Cauchy-Schwarz inequality, the second 
equality of Eq. (|A5[) implies that Vi — PijVj for all pos- 
sible pair with i < j. Here a complex number. 
The first equality of Eq. (|A5I) further implies every fiij 
must be positive. 

Due to the unitary invariance of the Hilbcrt-Schmidt 
norm, we have /3jj = 1 for all pairs (i,j) where i < 
j < k. When j = k + 1, since \\X\\ = \\vi\\ = 

/3i(fc+i)ii«fc+iii = pi(k+i)\\i:U A ^ XAi j\\ = 

based on Eq. (IA2[) . this implies that Puk+x) — 1 for all 
i < k + 1. 

Overall we have 

i 

UiXUl = u 2 xu\ = ... = u k xul =J2 a j xa ]- ( A6 ) 

Noticing that *(A) = XX, hence we get UiXLlj = XX 
for i G {1, 2, fc}, this completes the justification. 

Appendix B: PROOF OF LEMMA 3 

Proof. The adjoint operator of $ defined in Eq. © 
is given by 



&(p) = (l-q)tfpU + q U^l iP F xi U.(Bl) 

We shall first justify the following assertion: 

If #(A) = XX where |A| = 1, then &(X) = X*X. 

Theorem 1 implies that WXU — X*X and X — 
E l£Zw E ie{r , ;) Px,XPL- Note that <frt(X) = (1 - 
q)WXU + qJ2 x£ZN £ i£W} U^ xi XF xi U, this implies 
that &{X) = X*X. 

To prove statement 1, let Z G .E#(Ai) and Y 6 
•K*(A2). Then, by the assertion above, X\(Z,Y) = 
(&Z,Y) = (Z,<f>Y) = X 2 {Z,Y). Since A x ^ A 2 , it fol- 
lows that (Z,Y) = 0. Hence E^(X 1 ) _L £*(A 2 ). 

We proceed to justify statement 2. For an eigenvalue 
a with | a | < 1, we may assume, without loss of general- 
ity, that the generalized eigenvectors belonging to a are 
arranged in a sequence Yi, Y 2 , Yj a such that: 

(*-aI)Y r =y r _i, (B2) 

where, by definition, Yb = 0. It follows that (<fr — 
Q!l) r y r = 0. If Z G Ker($ - Ail), then, again by the 
assertion, X X {Z,Y X ) = {&Z,Y X ) = (Z,<f>Y x ) = a{Z,Y x ), 
which implies that {Z,Y X ) = 0. Similarly, X X (Z,Y 2 ) = 
(&Z, Y 2 ) = (Z, 3>Y 2 ) = (2", ii> + a(Z, Y 2 ), from which it 
follows that (Z, Y 2 ) — 0. Likewise, by the same reason- 
ing, applied repeatedly, we deduce that (Z, Y r ) = for 
any r. Thus E*(X X ) _L Span{Yi, F 2 , Yj a }. 

Appendix C: PROOF OF LEMMA 4 

Proof. Let A be an eigenvalue of 4? with |A| = 1. 
Picking X E Ker(*& — Al) such that X ^ 0, according 
to Theorem 1, we have £ x£Zjv E ieW} F^t/Af/tp^ = 
UXW which implies that UXW is a diagonal matrix. 
Since UXW = XX, then X must be diagonal. With- 
out loss of generality, it may be assumed that X = 

J2 X £i, N a xr\xr)(xr\ +Y,xez N a xl\ xl )( xl \- Note that U > 
as a linear operator, can be expressed by 

U = ^ [«n|(x + l)r)(xr| + M 2 i|(a; - 1)0 M 

-H*i2|(aj + l)r>(xZ| + u 22 \(x - l)l)(xl\\. (CI) 
Therefore it can readily deduced that 

UX = a xr un\{x + l)r)(xr\ 
+ a xr u 21 \(x - 1)Q (ar| 

+ ^2 a xiui 2 \(x + l)r) {xl\ 

+ a ^ u ^\( x - ^O^l; ( C2 ) 
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XU = ^2 a {x+1)r un\(x + l)r)(xr\ 
+ ^2 a(x-i)iu 2 i\(x -l)l)(xr\ 
+ ^2 a {x+1)r u 12 \(x + l)r)(xl\ 

+ X! a (x-l)l U 22\(x - 1)1) (xl\. (C3) 

Based on UX = \XU, we derive the following four 
identities. 

a xr u\i = Xa( x+ i) r un, a xr U2\ = Xa( x -i)iU2i- 

a x iui2 = Aa (2 , +1 ) r ui2, a 2;i U22 = Aa^-i);^- (C4) 

Here none of uii,ui2>W2i and U22 is zero, and x G Zjy. A 
little algebraic manipulation of these four identities im- 
plies that A 2 = 1 and A^ = 1. We arrive at the following 
two cases depending on the parity of AT. 

(1) If N is odd, then A = 1. By equations in (|C4[) . 
it is seen that X = 1&2N where k is constant, and so 
E&(1) = span{I 2 Ar}. 

(2) If N is even, then A=lorA = — 1. In the case 
when A = 1, X = kl2N where k is constant, and so 
E<&(1) — span{I 2 Ar}. In the case when A = —1, X = kl±i 
where k is constant, and Eq>(— 1) = span{I±i}. 

Appendix D: PROOF OF THEOREM 5 

Proof. In the case when N is odd, A= 1 is the only 
eigenvalue of 3? on the unit circle by Lemma 4, therefore 
the conclusion follows upon applying theorem 5 in (39j . 

In the case of N is even, by Lemma 4, we have E& (1) = 
Span{l2Ar} and E&(— 1) = Span{I±i}. In this context, 
we may choose for *B("H) a basis consisting of a basis 
Zi = ^^hN for £7*(1), a basis Z_ a = ^=I±i for 
Eq,{— 1), and together with an orthogonally complement 
basis consisting of generalized eigenvectors belonging to 
all other eigenvalues. In terms of such a basis we have 
p(0) = CiZx®C-xZ-i®W where W ± Z r for all r = -1, 1 



by Lemma 3. By simple linear algebra, it follows that 
c r = tr(Ztp(0)) for r = -1, 1. 

Relative to this basis, the Jordan canonical matrix rep- 
resentation of $ is given by: 

[#]=diag(l 1 -l 1 Ji,J 2 ,...,J h ) 1 (Dl) 

where J r is the Jordan block corresponding to the eigen- 
value a r of the magnitude strictly less than 1, r = 
l,2,...,h. 

Consider what becomes of the Jordan blocks of the 
powers [$]' as t y 00 . Since each of the Jordan blocks 
J r is an upper triangular matrix whose diagonal is popu- 
lated by a single eigenvalue of modulus strictly less than 
unity, it is a simple exercise in elementary algebra to show 
that lim t _ ) . 00 J* = O r (zero matrix of same size as J r ). 
Thus, if we define 

[*°°] (i) = diag (1, (-1)', Oi, 02, O h ) , (D2) 

then ||** - *°°(f)|| -> 0. 

Therefore ||*V(0) - ^h N ~ M) 4 ^(p(0), I±i>I±i || 
converges to zero; In particular, if (p(0), I±i) = 0, p(t) = 
**p(0) converges to j^hN- 



Appendix E: PROOF OF THEOREM 8 

Proof. We assume that N is odd. By The- 
orem 5, lim^oo p(t) = 5^l2jVj which implies that 
lim t _ i . 00 S(p(t)) = S(ojj^-2n) by Lemma 7. Since both 
trace w (•) and tracec (•) are continuous functions of the 
argument, it follows, by Theorem 5, that limt-j-oo pc(t) = 
trace w (^v^jv) and lim t ^.oo Pw{t) = trace c (^^iv)- 
These imply that lim^oo S(p c (t)) = ^(tracew (^f^iv)) 
and \imt-yoo S(pw(t)) = S*(tracec (^^jv)) by Lemma 
7. Because of S(^1 2 n) = 5(trace w (jfthN)) + 
5(trace c {^whN)), so lim^oo S (pc{t) : Pw(t)) = 0. 

In the case when N is even, we consider two sub-cases: 
(a) t is odd; and (b) t is even. In each of the two sub- 
cases, the assertion is proved in exactly the same way as 
in the case when N is odd. 
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